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Abstra ct 


This report represents the culmination of work to date in the areas of 
modeling and control of large space structures# Both theoretical develop- 
ments and the results of laboratory experiments are treated herein, as 
they apply to active attitude and vibration control, as well as static 
shape control, ^^odern control theory has been employed throughout as the 
method for obtaining estimation and control laws. 
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Section 1 

1.0 Introduction 

1.1 General Background 

As spacecraft becaite larger and more flexible, the equations needed 
to accurately model the static and dynamic behavior of these spacecraft 
became more complex. In general, higher order finite element models must 
be synthesized by the structural analysts in response to spacecraft tending 
away from lumped mass configurations and tending toward continuously dis- 
tributed configurations. More and more structural modes or degrees of 
freedom are passed on to the control analysts in an attempt to retain a 
faithful model. The use of these models directly for the purpose of on-board 
estimation and control of attitude, shape, and station may become unwieldy. 

On the other hand, serious well known stability problems may arise due to 
the use of oversimplified models. 

Spacecraft size also is not the only driver of model complexity. As 
more demands are made of control system performance, in terms of attitude 
accuracy, pointing accuracy, stability, shape accuracy, slew speed, etc., 
structures which heretofore have been modeled as rigid bodies must now 
include flexibility terms. This results from the fact that as more per- 
formance is required from a given space structure, improved knowledge of 
the structure Itself is vital, just as are more accurate sensors, actuators 
with better resolution, better computational accuracy and resolution, and 
a better understood disturbance environment. 

Prior to this work, very little information was available on the control 
of large distributed structures. That is, assuming that any sensor/actuator 
arrangement and type were available, the control logic for processing the 
sensor outputs to produce actuator commands had not been developed. The 
purpose of this work was to obtain a better understanding of the theory 
of the control of distributed parameter systems (DPS) through analysis and 
laboratory experimentation. Rather than attacking the analysis problem- 
by "building up" a general DPS from many discrete, Interconnected, lumped 
systems, and applying existing modern control theory, new optimal control 
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approaelies bosed on continuum models were developed for static am!? dynamic 
control. Using these now approaches, much Insight could be gained for the 
control of ger^'/al DPS modeled with finite elements. 

There arc many possible applications for spacecraft that can be modeled 
as DPS. Borne of these include antennas, solar arrays, platforms, and solar 
sails. Using antennas, a broad range of the electromagnetic spectrum con be 
covered. Wrap-rib, hoop-columu, and sunflower antenna designs may all 
be represented as DPS, and have applications in the microwave, radio, and 
K-ray bands. Solar arrays have applications ranging anywhere from power 
supplies for satellites to orbiting solar power stations for entire cities. 
Platform structures may provide conwion utilities Including power, thermal 
regulation, communications, and attitude control, for a number of different 
experiments. Finally solar sails may some day provide an economical method 
of developing low thrust over long time periods. 

In an effort to mimic many of the dynamic and control characteristics 
of DPv?, a laboratory facility has been assembled. The "spacecraft" consists 
of a hanging flexible beam. The beam is instrumented with position sensors 
and force actuators, A microprocessor directs the entire control sequence 
for shape control and vibration control. Laboratory demonstration and 
verification of various control system concepts will continue to be a vital 
part of the large space structure control effort. 

The work contained in this report on DPS can be divided into five 
major areas: 

1) Modeling 

2) Optimal Control 

3) Active Control 

4) Shape Control 

5) Experimental Verification 

The section on modeling consists of several of the more popular methods 
of mathematically representing physical systems for control system design 
purposes. The section on Optimal Control consists of a brief description 
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of some of the modern control approaches i Active Control refers to the 
control of systems governed by o set of dynamic equations of motion, This 
may include vibration suppression, attitude maneuvers, translation or 
station keeping maneuvers, etc. Shape Control means the control of an 
elastic structure to produce a given static shape. Finally, a section on 
Experimental Verification contains a description of the development and 
results of a facility designed to demonstrate and verify various control 
system design approaches. 
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Section 2 


2.0 Modeling ol: Large Space Structures 

2.1 Introduction 

The design of a control system for any physical system must certainly 
begin with some knowledge of the system itself. A mathematic.«l representation 
of the system, structure, or spacecraft, is known us a model, the modeling 
techniques un bo examined in this section arc partial differential equation 
models, finite element models, finite difference models, and modal models, 
bach type of model havS its own merits and shortcomings depending upon the 
particular application. In an effort to keep all tlie notation simple, 
undamped and non-rotating structures will be considered throughout the 
remainder of this work. Although tills type of model will not be universally 
applicable it will still represent a very large class of proposed large 
structures. 


2,2 Partial Differential Equation Models (PDE) 

The PDE is the most natural way of describing the behavior of a spatially 
continuous system, and, in fact, it will be the only continuous modeling 
approaci to be discussed in this work. PDE’s can be very useful for modeling 
simple systems, such as strings, beams, plates, diaphragms, shells, columns, 
etc. , but rapidly lose their usefulness for complicated structures. The 
render should be reminded t!iat there are many complicated structures that 
may be represented approximately using the simple models listed above. The 
primary value of the PDE model for control purposes is that a PDE is a very 
concise approach to handling a continuous model, i.e,, a model that retains 
an infinite number of degrees of freedom. Much insight into the control of 
general structures may be gained from the analysis of several "complete" 
models 0 


The general form for the PDE model considered here is given in (1) . 


p(x) y(x,t) 


y(x,t) -f f(x,t) 


B y (r,t) = 0 
y(x, t ) = y 



( 1 ) 
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The complete model consists of equations of motion, boundary conditions, 
and initial conditions. In other words, (1) can be stated as *'the mass 
times the acceleration of a y^int x in Q is equal to the applied forces." 
These forces are due to Internal forces and external forces. The dynamics 
are also governed by an appropriate set of boundary and initial conditions. 
The problem of static shape control may be incorporated into (1) by 
suppressing the time dependence. 

2,3 Finite Element Models (FB) 

The finite element modeling approach is a Lagrongian approach used 
to assemble the differential equations of motion for complicated structures 
from a set of simple elements, such as beams, plates, rods, point masses, 
etc. 


The basis of the finite element approach is that the kinetic and 
potential energies in a structure can be obtained as the sum of the energies 
of the individual elements, and that the energies In individual elements 
may be approximated using a (small) finite number of discrete coordinates, 
and some fixed continuous interpolation functions. The energy in a 
particular element results from a spatial Integration over the element, 
effectively allowing the energy to be written as a function of these 
variable discrete coordinates, and of some fixed constants that result 
from the integration process. Applying Lagrange's approach for deriving 
the equations of motion results in 

Mx + Kx * F 

Attitude control, stationkeeping, and shape control of many structures 
can be incorporated in the format in (2). 

The advantage of the finite element method of modeling is that models 
for arbitrary structures may be synthesized from much simpler component 
pieces. The mating of various components, as well as the boundary conditions 
of the overall structure are satisfied "automatically." Furthermore, a 
great deal of software has been developed for assembling high order models, 


and for performing the subsequent analysis. One disadvantage of finite 
element modeling Is that the model itself may be very high order, 
(perhaps hundreds), and that useable results can only be obtained 
after an eigensystem analysis of the high order model. 

2.4 Finite Difference Models (FD) 


A finite difference model is a direct approach to deriving equations of 
motion, as oppose i to the variational approach analyzed with FK. The 
finite difference method is the result of approximating a differential 
operator directly by means of finite differences. In general, therefore, 
the higher the ord^r of the operator, the more coupled the dynamics of a given 
point is with adjacent points. General boundary conditions can only be 
enforced by painstakingly obtaining the FD equation that is also consistent 
with thi; desired boundary value. It is not possible using FD techniques 
to as jcmble complex models using simple components, as is the case with 
F12. FD modeling for large structures lias only very limited usefulness. 

The general form for a finite difference model is the same form as (2), 
except that the mass matrix with FD is diagonal. 


2,5 Modal Models 


A modal model of a structure is usually the output of another modeling 
process, generally PDF, FE, or FD. Modal models yield a decoupled set of 
second order differential equations describing the natural or unforced 
behavior of the structure, in terms of the modal amplitudes. Because of the 
fact that the natural modes of a system are easily observed, much Insight 
can be gained by study of a modal model. Under the assumption that any 
motion of the system can be represented as some time varying linear 
combination of the modes, knowledge of the mode shapes and their time varying 
amplitude is sufficient to completely describe arbitrary behavior. The 
Inhomogeneous modal model of the time varying amplitudes is usually written 


q + 


‘n 


<ji^(x) f (x,t) dx 


(3) 


where the eigensystem analysis is either performed analytically in the case 
of PDE, or numerically (using a wide variety of available software) in the 
case of PE or FD models. 
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For tho aforementioned reasons, It would appear that modal models are 
ideal for controls applications, being natural and simple in structure, 

In fact, modal models had enjoyed widespread usage up to and including 
their application to large structure control. The particular problem 
encountered with the use of modal models is that in theory, an infinite 
number of modes are required to model a continuous structure, whereas in 
practice, it is only possible to work with some finite number. Truncation 
of the series approximation of the model may result in system instability, 
inasmuch as these modes are still forced by outside forces, but the 
resulting motion was previously assumed to be unimportant. 
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Section 3 


3.0 Optimal Control 

3.1 Design Procedure 

The primary reason for adding a control system to a system or a structure 
1b to change the characteristics of that system, whether it be in response 
to commands, or in response to disturbances. Kvidently, there arc changes 
in the system that can be made which makes its characteristics closer to 
some desired performance, and there are other changes that make the perfor- 
mance loss desirable. One method of specifying a measure of performance in 
the control system design process is through control system design methods 
using optimal control. 

The optimal control design process involves specification of a scalar 
performance measure or cost function that reflects the relu/ivc Importance 
of a variety of different factors. Tiie choice of feedback control law that 
minimizes the cost function is then deemed to be "the optimal control." 

The form of the cost function used throughout this work is a quadratic 
performance index. That is to say, a quadratically increasing cost is 
incurred for both (1) using additional control effort, and (2) desired outputs 
from the system not matching actual outputs. 

Quadratic cost functions yield easily implemented linear state estimators 
and controllers. Furthermore, an entire body of information exists for the 
theory, performance, and implementation of this type of controller. 

3.2 Sensors and Actuators 

The purpose of this work is to address the theory of the control of 
large structures. The two major reasons for the need of a new theory of control 
is that future large space structures differ from present day structures 
in two very significant aspects. First, the control of large structures 
involves the control of systems governed by continuum or infinite dimensional 
models, and second, even a simplified model of the structures is likely to 
contain gross errors due to the inability of conducting meaningful ground 
tests of these large structures. 
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The problem of deelgnlng new sensors nnd ncuviators, ov, of selecting 
them from OKistlng senaor/actviator technology Is not addressed in this 
work. When needed, it will be assumed tluit ideal position sensors, angle 
sensors, rate sensors, uorquers, etc. will be available for the control 
system design process. 

3.3 On-Board implementation 

As with the sensors and actuators, very little attention will bo paid 
to the on-board implementation of the control system. It is implicitly 
assumed, however, that an on-board dltlgal computer will be available for 
implementing the control law, if necessary. In a qualitative sense, it 
will be recognised that there are limitations to the computational capability 
of the on-board processor, whether the processor is physically in the form of 
a single, lumped computer, or if it consists of a distributed matrix of 
smaller processors. 
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Section 4 


4.0 Active Control 

4.1 Introduction 

Tlii« section contains a derivation of the necessary conditions for an 

optimal local control law for a general Byatc'm in state variable format. 

A local control law is one in which only local state Information is used 

to synthesize the control law for each actuator (even If additional state 

information is available) and hence, is a different problem from that of 
(2 3) 

output feedback. * A practical application of this idea occurs in 
designing control systems for large flexible space structures, where many 
sensor outputs may be available for feedback, yet it Is not practical to do 
so due to the large spatial distances Involved. 

by representing the structure, in physical coordinates (the initial out- 
put of the PDU,l’l) or FK analysis) rather than modal coordinates, the numerical 
solution of the optimal full state and local feedback problems may be 
simplified. Numerical examples of control law designs for a simple two 
mass model, for a free- free flexible beam, and for a string in tension 
are given. Similar examples using different design approaches may be found 
in References 4 and 5. 

The counterpart to the local control law is a local estimation scheme. 
Local estimation processes sensor Information using a dynamic model of the 
system to obtain estimates of only the nearby components of the state 
vector, not the entire state vector. These compocents of the estimated 
state vector should be precisely chose components required by the local 
control law* Combining the local state estimator with the locAl control 
law results in a local controller wUlcli significantly reduces the amount 
of on-board computation, and allows the computations to be performed in 
a distributed or parallel manner. 

4.2 Necessary Conditions for Optimal Local Control 

The derivation of the necessary conditions for the optimal local control 
begins with a system in state variable format and a quadratic performance 
index to be minimized as is shown in (4) . 
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X « Fx Gu + w 


Min J ■ 
u 







(x*^ Ax + u*^ Bu) dt 


given 




Assuining that the desired final solution is of the form u •• Cx, the 
substitution into (4) is made yielding the deterministic equivalent 


(F + GO) X 
/■t 


x(t^) 


given 


(5) 


Min J 
C 


L/2 


^ BC) X dt 


0 


At this point, the constraint that only local states be fed back can be 
enforced by requiring that certain components of C be identically zero and 
that the minimization in (5) be carried out wltli respect to the remaining, 
non-zero components of G. Equivalently, the constraints may be adjoined 
to the Hamiltonian to yield 


H - 1/2 (A + c'^ BG) X + X^' (F + GC) x + <b) 

where 

^ij “ ^ ^ij ^ 

^ij ° ^ij “ ^ 


and where summation over the repeated Indices is implied. The optimality 
condition becomes 


T T T 

BC xx^' + G X X + u ■ 0 


(7) 


where the picked to make the corresponding constrained 

equal to zero. 


The usual sweep solution obtained by letting 

X(t) » S(t) x(t), C * - B"V S(t) (8) 

no longer works in general since the optimality condition in (7) results 
in C being a function of x and t, i.e. the minimization can no longer be 
performed independently of the initial condition. However, the appearance 
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T T 

of Che terms xx and Xx above ouggeats that a linear statistically optimal 
control may exist. The stochastic analog of (5) in terms of the covariance 
of the state is^^^ 


X « (F + GC) X + X (P + GO * + Q 

ft 




given 


( 9 ) 


Min J « trace 
C 


^ 1/2 (AX + BC X) dt 


0 


.Ts 


whore deeignated are zero and X » E(xx ). The accomplishment In the 

preceding step is to average the performance index over a range of possible 

initial conditions. Rather than considering all initial conditions to be 
C7) 

equally likely, a more realistic range of possible initial states can be 
obtained using the state covariance matrix. The adjoint matrix equation 
and optimality condition ares 


-A » A (F + GC) + (F + GC)'^ A + A + c''-' BC A(t^) * 0 
BC A + 0^' h X + y - 0 . 


( 10 ) 


Although an exact solution to (10) is possible, an approximate solution may 
be easily obtained by expanding the equations in (10) to first order in 
y about the optimal solution for y * 0. This yields 

C(t) a - b"^- [g'^’ A(t) + y X"^(t)] (11) 


where A and X are the solutions of the unconstrained optimal control problem 

(y » 0), and y is picked to zero the corresponding components of C. For 

the unconstrained problem, y « 0, the result in (11) reduces to the familiar 
T 

result, C*'-BG A. Equation (ll) has a nice physical interpretation. To 
first order in y , if not all the states can be fed hack, those states that 
are available should be fed back with a correction to the feedback gains 
based on the correlation between those states fed back, and the remaining 
states. (As a practical note, the inverse in (11) need never be computed. 

In fact only "a few" elements of X need to be manipulated.) Furthermore, 
since the solution is expanded about the optimal solution, it can be shown 
that “ 0, i.e. this concept of local control does not severely affect 

performance to first order in y. The algorithm for solving the local 
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control problem can be outlined as followi. First, »olve the full stats 
optimal control problem, next, apply the local control correction appearini; 

In (11). 

It should be noted that as in the full state feedback case, for F, 0, 

A, B, and Q all constant, it is possible that a steady state solution for 
C, A, and X may be obtained ns tj. - t^^ ► «. However, as opposed to the case 
of full state feedback, stability of the closed loop system is not guaranteed 
wiien using local control gains. The eigenvalues of F + tiC using the local 
control gain must be determined to verify suability of the closed loop 
system. 

Thus far, the problem of solving the full state optimal control problem 
(step one above) for high order systems has been avoided. However, the 
solution procedure required usually aasuraes that a low oi’der model of the 
high oi’der system Is available to make the problem tractable. Traditionally, 
the structural analyst supplies as many modes as the control system Ueslgner 
wishes. Regardless of where the actual truncation occurs, or by what method, 
the control system designer begins with a deficient model. This may result 
in closed loop instabilities. 


The alternative to working with the truncated system model ivS working 
with a full order finite element model, as it is generated by the structural 
analyst. The next section offers some hope that the analysis of these 
high order systems, particularly structural systems, may be feasible, 

4.3 Finite Element Structural Models Control Design 

In an attempt to develop control system design techniques for high order 
systems, and to alleviate the problem of truncated modes, it is worthwhile 
to examine the structural equations of motion. Attitude control, station- 
keeping, and figure control of many structures can be represented by the 
following matrix equation. 

Mx + Kx - Gu (12) 

Matrix handedness of M and K is a direct result of the finite element modeling. 
Since control inputs from a given actuator are applied at a single station 
on the structure, nearly all of the elements of G are zero. 
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Equntion (12) can be placed in «wdal form 


" 2 T 

q 4* 0 q » Gu 

by normalising che eigenveetor matrin, , ho that 


( 13 ) 


X * ^ q 

M X (lA) 

<!»* K « • U*’ 


where U 1h a diagonal matrix of modal frequencies. Furthermore, computer 
programB like EIGSOl/^^ and can solve Che open loop eigenvalue 

problem (u » 0 in (12)) very efficiently by taking full advantage of the 
matrix aparsity in both storage and computation. 


Consider now the problem of designing a control system for the system 
in (13). Selection of a control law can be based on the minimization of a 
quadratic performance index similar to that appearing in (4) . 


J " 1/2 



(x^ Ax 


■f u Bu) dt 


(15) 


Regardless of where or if truncation of the modal system occurs, the open 

2 

loop system dynamics matrix, Sf , is a diagonal matrix and hence, simple to 
manipulate from a computational point of view. However, the corresponding 
control distribution matrix (j) G s not a sparse matrix, and so the 
Hamiltonian system for the corresponding optimal control problems has the 
following form 



The shaded areas in (16) represent non-zero matrix entries. Because little 
useful matrix structure remains in (16), the eigensystem analysis required 
for the solution of the optimal control problem can not be performed 
efficiently. Consider Instead the control of the original dynamic system 
in (12). At this point the concept of local control emerges naturally. 

The dynamics of the flexible structure are characterized locally. This is 



the reason that a good dynamic model of the flexible structure can be 
obtained using tightly banded matrices. Furthermore, actuators produce 
effects locally, and sensors measure local behavior. It therefore seems 
plausible that a good controller may be possible using only local state 
information. 


The optimal control problem may be formulated as follows. 


■’ t 

Mx + Kx » Gu x(t^), x(t^) given 

ft. (17) 

Min vJ » 1/2 (x^Ax + ii Bu) dt 
u t 

J O 

By adjoining the constraints in (17) to the performance index with Lagrange 
multipliers, X, and integrating by parts, two times, the closed loop system 
dynamics and the corresponding matrix structure are given by 



x(tjj), x(t^) given 

x(tp - X(tf) - 0 


Equation (18) is the Hamiltonian system corresponding to the optimal control 
problem in (17). The important feature to notice is that nearly all of 
the original system’s matrix structure is preserved, and that the eigen- 
system analysis that must be performed in (18) to obtain the optimal 
control differs from that which must be performed in (17) to analyze the 
open loop system dynamics, by only two ’’stripes" off the principal diagonal. 
As such, from a computational point of view, both storage and computation 
time can be reduced by fully exploiting the high degree of matrix sparsity 
in an eigensystem analysis or through efficient matrix perturbation 
techniques. 


4.4 Examples of Local Control Systems 

This section contains two examples of the steady state optimal local 
control concept and comparisons of the performance of the optimal local 
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control law with the performance of control laws designed using various 
other approaches. 

Example I - Two Mass Model 



min E(J) 
Fig. 1 


Kx2 

J • 1/2 j + 4u^ 

^ 0 

Two Mass Model System 


dt 



Figure 1 consists of two unit masses connected by a spring with stiffness, 
one half. The open loop system has a single rigid body translation mode 
and a single vibration mode. It is desired to control the position of mass 
1 (x^) with a control input u, in the presence of the disturbances w. 


Three approaches for designing control systems are examined; 

a) Full state optimal control. All four states are available for 
feedback (x^, x^, x^, • 


b) Modal control. The control system is designed on the basis of a 
rigid body model. 

c) Local control. The feedback gains on X 2 and X 2 are constrained 
to be zero. 


For the given performance index, the results are most easily summarized 
in figure 2 and Table 1. 

Table 1. Control System Design Results .. 


Controller 

eigenvalues 
Mode 1 Mode 2 

*1 

gains 

*1 *2 

h 

Performanci 
Index J 

Open Loop 

6.0 

± 

0 

0 

0 

0 

m 

Full State 

-.310 + .404j 

-.128 ± .974j 

-.384 

-.8/6 

-.116 

-.453 

3.792 

Rigid Body 

-.357 + .619j 

-.355 + .606j 

-.500 

-1.414 

0 

0 

4.949 

Local 

-.144 + .299j 

-.186 + .984j 

-.221 

— .661 

0 

0 

4.393 
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X OPEN LOOP 
• OPTIMAL FULL STATE 
A MODAL, DESIGN POINT 
□ MODAL, ACTUAL LOCATION 
L LOCAL 




X- 


1 


Figure 2. Control System Design Results. 



The closed loop eijjenvalues, feedback gains, and performance index 
are given for each control system design approach. Of course*, the full 
state optimal control law ’’performs" the best. The desired root locations 
for the control system dcvsigned using the truncated dynamic model are shown 
in Figure 2 with triangles. Due to the presence of the truncated system 
dynamics (in this case tlie vibration mode) , the four closed loop roots 
actually end up some distance away from the rigid body design point. 

The local control system is designed using the full systems dynamics, but 
with the constraint of partial state feedback. This closed loop system is 
’’closer" to the optimal full state control law than is the rigid body 
control law in terms of both the system performance and the final root 
locations, and only uses feedback of and x^^. It should be noted that 
as the cost of control decreases (B decreases), the relative merits of 
the local control law over the rigid body control law become even more 
apparent. 

Example II - Free-Free Flexible Beam 


The partial differential equation of motion for a free-free flexible 
beam with constant properties per unit length is given below (fig. 3). 



Position 


3^y . 
P — V + 

El ^ - u (x,t) 

^(0). 

.2 

■H(L) 

at 

3x 

3x^ 




^(0) . 

■^(L) 



3x^ 

3x^ 


m 0 


- 0 


Fig. 3. Free-Free Flexible Beam. 

A finite difference discretization of this beam can be obtained by 
choosing a state vector composed of deflections and deflection rates at 
ten stations along the beam (fig. 4). 



Fig. 4. Dlicretlzed B<*«m Model 

Furthermore, It is assumed that control forces u^ can be applied at the 
designated stations. Penalizing the beam displacements at each of the ten 
stations results in a control law which performs stationkeeping, attitude 
control, and shape control. The following matrices were chosen for the 
quadratic performance index 

A » B - 0.01 • 

As before 

^ “ ^ 10 * 

A full state optimal control law and a local control law were designed 
for the flexible beam using p « El =* 1, L « 9. Due to the symmetrical 
placement of tlie actuators, it is sufficient to present the feedback gains 
for synthesizing u^ and u^. In each of the accompanying figures (figures 5,6,7, 
and 8) solid lines are the full state feedback gains as a function of station 
location and the broken lines are the local feedback gains obtained under 
the constraint that only states which are immediately adjacent to the 
actuator are allowed to be fed back. 

There are several Important features of these results to recognize. 

1) The full state feedback control law makes very little use of 
"distant" state information. This feature is not apparent from the modal 
approach where the feedback gains corresponding to the various modes may 
be roughly the same magnitude. Evidently the modal feedback effects tend 
to accumulate at the actuators and cancel far from the actuator. 
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FEEDBACK GAIN (FORCE/DISPLACEMENT) 



Figure 5. Full State and Local Position Feedback to Actuator 1 


20 



Figure 6. Full State and Local Position Rate Feedback to Actuator 1. 
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Figure 7. 


Full State and Local Position Feedback to Actuator 2 


FEEDBACK GAIN 


iiGCNOt 

tOCAL 

WUHATE 


STATION ALONG BEAM 





Figure 8. Full State and Local Position Rate Feedback to Actuator 2 
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2) The local control law Is equivalent to active springs and dasUpots, 
to provide stntionkceping» attitude control, and figure control. This 
result may be useful for control system design for general flexible 
structures. 

3) In this example, the difference between the performance index 
using the local control in place of the full state optimal control is 
less than 1%. In general, the use of local control guarantees no first 
order change in the performance index. 

4) From a computational point of view, the following computer CPU 
times were required by a UNIVAC 1108 to obtain the control laws for this 
20-state system. 

a) OPTSYS,^^ 18 seconds. The QR algorithm is applied to the 40th 
order Hamiltonian system to extract closed loop eigenvalues and 
eigenvectors. 

b) Direct Integration, 15 seconds. Direct integration of the matrix 
Ricatti equation exploiting all matrix sparsity was employed to obtain 
steady state gains. 

4.5 Local Control Based on PDU Models 

One possible method of circumventing the problems resulting from the 
increased number of modes to be controlled is by avoiding modal models 
completely. Working directly wit.i partial differential equation (PDE) 
models of continuous spacecraft is a viable alternative. In fact, some 
Structures, such as solar panels, large antennas, and astromasts may be 
more easily modeled with PDE than with FE, particularly when many modes 
are required. 

The motivation behind the "local" approach to be discussed here is 
that the PDE for a continuous spacecraft describes the acceleration of 
eacli physical point on the spacecraft in terms of differential operators, 
i.e. in terms of the behavior of the spacecraft within a local neighborhood 
of each point. Since sensors measure and disturbances affect local variables, 
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a reasonable job of state estimation may be accomplished with local state 
estimation. The entire design procedure would avoid entirely the problems 
associated with high order (possibly truncated) modal models. If it la, 

In fact, true that a local control law minimises a quadratic performance 
Index for distributed parameter system, then the underlying explanation for 
this certainly must be derivable from the PDE control formulation. 

The control analog of this problem has already been addressed In the 
previous section. It has been shown that the optimal control law for a 
free-free flexible beam very closely resembles a local controller. This 
means that if a modal controller for this system Is designed, the feedback 
law may actually represent a more simply expressed control law were it 
expressed in physical coordinates, 

( 11 ) 

Breakwell' has also obtained some results for the optimal control 
of a continuous flexible beam using PDE, in which a symmetric root locus 
approach is used to determine the optimal closed loop root positions. 
Following this step, tlie continuous control gains needed to move the roots 
to their optimal closed loop locations remain undetermined. That is, no 
direct procedure is available for directly determining the optimal feedback 
gains. 

The purpose of this section is to present a PDE estimation and control 
problem formulation that contains sufficient generality to encompass a wide 
variety of continuous control system design problems within a single 
analytical framework. By drawing on a simplified example, a string in 
tension, some insight to the control of continuous structures is obtained, 
and some generalizations for their control can be made. Although the 
estimation problem is discussed here, the results can be extended easily 
to include the control problem too. A procedure for the direct determination 
of the continuous optimal feedback gain is given, following the format 
used in the control of systems governed by ordinary differential equations. 

Consider a general partial differential equation of motion in state 
variable format. This equation represents the dynamics of the state vector 
in some spatial domain, with boundary conditions prescribed on the boundary 
of this domain. The initial state of the system is given. It is desired 


to obtain an eetimate of the «tate of tho iyitem aa tlm® evolvos, from 
noiay meoBui'eraonts in the presence of cliBturbancQS. Th« mathematical model 
if given below. 


y(x,t) - P(D) y(x,t) 4* r(x,t)w(x,t) 
y(x,t^^) given 

b(n) yd’it) given 

z(x,t) *• U(D) y(x,t) + v(x,t) 



(19) 


F(D) , n(D), nml B(D) are Uncar differential operators in the spatial 
variable. Were it not for the preaencc of the differential operators in 
(19), the mathematical formulation would be Identical to the ordinary 
dlfferentinl equation state variable format* (It will theroforc be nota~ 
tionaXly convenient to omit future appearances of all superfluous symbols. 

It should be remembered that the model actually analyzed is written in all 
of its detail only in (19)). In order to obtain a state estimate from the 
measurements, a performance criterion ie selected which weights the relative 
importance of the process and measurement disturbances. Minimization of tills 
porformanee index is one way of obtaining the optimal state estimate, 


J - X/2 


ft. 


(x^^Q ^w + v*'!! ^v) dx dt 


( 20 ) 


0 


Adjoining the constraint dynamics in (19) to the performance index 
irv (20) gives 


ft. 


[ -f (z-Hv) ^'ir-^(z-Hv) 

0 ^ 


,T„-1. 


( 21 ) 


+ X (-y + Fy + rw)} dx dt 
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The neceiiory condition for the minimisation of the performance index 
ia that the flrat variation of J vanlahea for arbitrary admissible variations. 
This variation yields.* 


0 « W 

ft. 


n 


(wV^ + x'^’r) 5w dx dt 


{-6yVR*’^(z-Hy)-X^’ay 


+ X F ay) dx dt 


( 22 ) 


Since aw in (22) is arbitrary, the optimality condition ia obtained. 

w - - Q r’^X (23) 

Two steps are now required to isolate the ay multiplier in (22), 

T * 

The first is to integrate the X 6y term by parts to obtain the final 

T 

conditions on the adjoint variable, X (x,tr) » 0, and also an integrated 
• T ^ 

term, -X ay. The second step requires the isolation of the 5y terms from 

T 

their operators F and U . A short digression is required. 

Let (•»•) define un inner product. 


If T is some operator then its adjoint operator T* is defined as 


(u,Tv) « (T* u,v) 

Using the star notation to denote adjoint operators, the equations for 
the adjoint variable becomes 

X + F)*fX - H’^'^ir^(z-Hy) - 0 (24) 

Boundary conditions on X will result from computation of the adjoint 
operator. 

The closed loop dynamics of the estimator is now written below. 


1 

m 


■» 

r* •- 


F -rqr^’ 


y 

i 0 

m 




+ ■ 


-F* 


X 

' T -1 
: tr*R H 


* 



t 


y <t„) 

given 

By « 0 

(25) 

X (tf) 

« 0 

B*X » 0 
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( 26 ) 


As uguoX, It is dasirablt to writs ths sstimator in ths form 

}. 

y * Fy + K (K-lly) 

Hsrs K inoy bs a cllffsrsntiol opsrato? asti»n«tor gfllrit Toward this and, 
let 

A 

y ■ y « PX 

where P is also an operator. It is now straightforward to show 
y « Fy + K (t“Hy) 

p «• PF* -p FP + rqr^' - hp p( 0) « o C27) 

T -1 
K « piri*fR 

where products of operntore of course represent composite oparationr. 

The problem of determining the opilW'^l estimator gains for a spatially 
continuous system has now been reduced to familiar format# The only 
differences are that the estimator gains are now in operator form, and 
that the RIcatci equation for the estimator gain is in terras of a differential 
operj,tor, P. That such a solution exists and is unique is not pursued 
hero.^^^^ Rather, a detailed example problem is presented in the following 
section, 

A. 6 Flexible String Example 

Consider the dynamics of the string in tension shown below. 


The partial 


2 

iJL 

at^ 



differential equation of motion is given 


X Position 
by 


-p w y(0,t) " y(l,t) “ 0 
3x 

y(x,0) » 0 


( 28 ) 
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Assunia a continuous position measurement, z, is avollablo, and that 
it is desired to find the steady state estimator that will minimize the 

following performance index, 

,00 , 


1/2 




<w^q“^w + dx dt 


(29) 


From (28) and (29), the equivalent state variable formulation is obtained as 


F * 

Letting 

Cu,v) 


0 1 
0 

0 ^0 


i H « [1, 0] . 


(30) 


1 : 


uv dx dt 


the adjoint operator to F is given by 


(13) 






0 


(31) 


Using the results in (27) and letting 


2i 


P3 


V! 


it can be shown that the steady state optimal estimator can be obtained 
from solution of 


P3 + Pj - P^R-lp^ . 0 
P, + P.D^ - P^r“^P,; » 0 

Si X / ^ 32 ) 

D^P^ + P^ - ■= 0 

D^?2 + 230^ + Q - P3R~^P2 = 0 

Again, it is pointed out that products above denote composite opera- 
tions. Even so, the feedback gains can be written symbolically as 


K = 


r '’P" 

I P3/R 



(33) 
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where the square root of an operator is to be Interpreted as the operator 
whose composite operation with itself yields the desired operator. In 
each case, results obtained using the above analysis will be compared 
wit]) mtmericnl results using a discretized 20 state representation of the 
string dynamics, with the program OPTSYS.^^^' 

In the following examples, the estimator gains determined using the 

4 

differential operator approach will depend on the relative norms of D 

4 

and Q/U in (33) . Obviously the norm of D defined as 

" (y.o^y)/(y,y) - (D^y,o^y)/(y.y) 

can be made arbitrarily large. At this time some notion of the modal 
concept does prove to be useful (although it is still not required). In 
general, highe.r frequency modes in the time domain result in higher frequency 

4 

mode, shapes in the spatial domain. The norm of D on these eigenfunctions 
will then serve as a basis for which modes are low frequency, and which are 
high frequency. What will eventually be seen is that there will be 
dlstlnet estimation approaclies for the high frequency modes versus the low 
frequency modes. 

Case I Q »» lO'^, R » 1 

For Q/R large, or alternatively, for estimation of the component motion due 
to low frequency behavior, an expansion of (33) gives 



Regardless of the dynamic system involved, for Q/R large, the estimator 
gains as a function of position are approximately constants. Since they 
are not "operators”, they will involve feedback of only local state 
Information. In fact, Q/R large is precisely the condition needed to 
guarantee the optimality of local estimation (or local control in the 
control problem) . The estimator for this problem will then be 
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iL 

dt: 


(35) 


y 

X 

y 


0 




*^1 


(«-y) 


Notice that the estimate crrot tends to zero for all modes, not just for 
the low frequency modes. That Is to say, truncation of the differential 
operator to include only the leading term does not result in any instabilities. 
A compariBon of the estimator gains obtained using the operator approach 
versus the discrete approximation is shown below (fig. 9). 


The real power of the operator approach can be seen in the next 
example where the sensor accuracy spatially varies. 

Case n Q » lo'^, R « /37x 


In this case, no measurement is available at x*=0, and the accuracy of 
the measurement Increases as x increases. The norm of Q/R is identical 
to the previous case, so the same expansion of the differential operator 
applies. In this case, tl\e first order operator approximations to the 
optimal estimator gains arc 



(36) 


Again, the feedback gains obtained via tlie operator approach are compared 
to those obtained via the discrete, approximation (fig. 10). 

By allowing R to be a function of x, it Is possible to approximate the 
case of spatially discrete sensors by choosing appropriate continuous 
approximation functions for R(x) . 

Case III Q = lO'*'^, R ^ 1 

For Q/h Small, or alternatively, for estimation of the component of 
the motion due to high frequency behavior, a different expansion of (33) 
gives : 
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Figure 10. Operate-^ vs. Discretized Control System Design Results 
(variable weight) 
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( 37 ) 


l-j/K 


2 ^ 6 
21) 8ir 


Pj^/R 


M /■■ 


As would be expected, for Q/R * 0, the steady state estimator gains 
are zero. For Q/R not equal to zero, tlie inverse powers of D must be 
interpreted as Inverse operators. For the case of the flexible string, 
the inverse operator to I)^ that satisfies the boundary conditions is 
simply the integral operator whose kernel is the Green's function for 
the string, i.c. 


K f « g.(x,f,) f(r,) d? 

i r ■' 0 ^ 


(38) 


where 




|S(x*-l) 

|x(5"l) 


X > r, 

X < f. 


Not" 


that -1/D^ 


has a positive norm. 


The analytical expression for the /^/D operator is sc lewhat more 
complex, and somewhat more interesting. 



g (x,0 HO dC 

0 


where 


JL sin ir(x+g) 

TT sin r(x-5) 


(39) 


It should be noted that this operator satisfies the required boundary 
conditions at 0 and 1, but has the peculiarity of being infinite at x=5. 
At first, it might not seem reasonable that a particular point along the 
string could have an infinite weight contribution to the control law. In 
fact, this is not the cose. The kernel does have an infinite value at 
x<=5, but the fact that this is an integrable singularity, and that the 
kernel does appear under an integral, puts this feature in the proper 
perspective. It also explains why the gains computed using a discretized 
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system dynamics are not infinite valued at but instead ore an average 
value over the discretization interval (fig. 11). 

The Green's function and its square root yield control and estimator 
gains and appear naturally in the ?DE formulation for controlling (damping) 
high frequency modes. These concepts can be extended to Include generol 
models based on EE by using the elgensystem expansion of the Green's 
function. 


The Green's function can be expanded as 

g (x,e) - I 1 

^ i-1 


and 


» tj). (5) 

“ I —71 

^ 1«1 


(40) 


where the t{}'s and X's are the normalized eigenfunctions and eigenvalues 
of the operator. For the case of an FE control design, the elgensystem 
analysis is usually performed as part of the structural analysis. The 
fact that these results can be directly used for control system design, 
as well, is a valuable nev^; result. 

The results obtained in this section are generally applicable to a 
wide range of PDE models, not just the string in tension. The general 
result can be summed up as follows. 

Local estimation and control Is optimal for the higli gain, low frequency 
modes, where performance is required; and estimation and control based on 
Che Green's function and Its square root is optimal for the low gain, high 
frequency modes, where damping augmentation is required. 
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Section 5 


5.0 Sliape Control 

5.1 Introduction 

Large, lightweight, flexible, space structures exhibit dynamic shape 
variations greater than those of any previous spacecraft. Whereas the 
shape varies continuously with changes in space and time. It must be 
estimated (controlled) by sensors (actuators) placed at discrete points along 
the structure. 

The. mixing of continuous and discrete mathematics can present or ease 
difficulty in the solution of control and estimation problems, depending 
on the method by which it is approached. The use of Green's functions to 
convert boundary value problems into integral equations provides a 
convenient treatment of these problems, 

The Green's Function provides a solution to a nonhomogeneous ordinary 
or partial linear differential equation in terms of an integral operator 
which acts on the forcing function (the nonhomogeneous term). Comparison 
of solutions for different forcing functions becomes relatively easy. 
Furthermore, the expression of th? solution in terms of the Green's 
function is especially convenient in the case of the spatially discrote 
functions found in large space s true tun- control and estimation problems, 
since some integrals become finite sums. 

The use of an integral operator rather than a differential one 
possesses additional advantages; 

(1) The expression of a solution as an integral equation automatically 
incorporv-ites the boundary conditions, which must be handled separately 

if the problem is stated as a differential aquation. 

(2) The integral operator is usually bounded and often completely 
continuous, whereas differential operators are unbounded. Tlius results 
concerning eigenfunction expansions, solutions of nonhomogeneous equations, 
etc. ate more easily obtained. 
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(3) Numerical approximations and variational techniques which Include 
several other methods of solving problems with constraints are more easily 
applied to integral rather than diffei*ential equations. 

5.2 The General Boundary Value Problem and Green’s Function^^^’^^^ 

We first define a general boundary value problem and discuss the use 
of the Green’s funation in its solution. We then apply the technique to 
the solution of a goner^X control problem and a general estimation problem 
for large space structures. 

Consider u surface whieli occupies a simply connected region i] and la 
bounded by the curve F. 



r 

Assume the surface is acted on at each polntPciJby a force f(P) and that the 
static deformation u(P) of the surface satisfies the partial differential 
equation 

Uj - f (Al) 

where L is a linear partial differential operator, and also satisfies 
appropriate boundary conditions 

B^(u) « 0, 1 1 i 1 N> for Per . (42) 

Assume the boundary conditions (42) are such that the operator L is self- 
adjoint, That is 

(Lu,v) « (u,Lv) (43) 
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for any pair of functions (u»v) In an appropriate class which satisfy 

the boundary conditions (42), The Inner product (u,v) is defined to be the 

Integral 

(u,v) » f u(Q) v(Q) dQ . (44) 

hi 

We also assume, for convenience, that the homogeneous system 

Lv “ 0, B^(v) « 0 , 1 < 1 1 N (45) 

has only the trivial solution. This is equivalent to the assumption that 
the system has no rigid body inodes. 

Thus, the general boundary value problem is the system (41-42) 
together with the assumptions (43-45), 

The Green's func...lon g(P,C)) associated with the boundary value 
problem (41-42) satis flOvS 


hgdSQ) - 5(P-Q) (46a) 

B|(g) * 0 , i 5 i i N , for Per , (46b) 


It 

Q. 

to 


represents the response at the point P to a unit impulsive force at 
Since L is self-adjoint (u,].p^ “ (Lu,g), which provides the solution 
the, boundary value problem (41-42) : 
r t 


u(P) 


u(Q) S(P-Q) dQ - 


g(P,Q) f(Q) dQ . 


(47) 




Hi 


Remark 1. The Green's function is the kernel of the compact Integral 
operator K such that 


U 


t 

g(P,Q) f(Q) dQ . 


(48) 


K is clearly the inverse of the operator L, where defined on the range of 
1,, since KLu s* Kf “ u and LKf =* Lu “ f. 


Remark 2, The solution of (46a) is called a fundamental solution. The 
equation (46a) is satisfied in a distributional rather than a pointwise 
sense. That is 

(hg,<ji) « (g,L*4>) “ (^9) 
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for iiU test functions (A fimation Is on infinitely cliff erentlabl© 
function defined on a domain which has compact support,) 

There are some additional requirements for a solution in the ease 
that (*i5) has hon-trivlal solutions, that is, the physical system possesses 
rii^id body modes. • In this case, (^il) has no solution unless the "consistency 
condition" 

(f, v) « 0 

is satisfied for every v(Q) which is a solution of (45). 

The consistency condition becomes reasonable when we consider that 
seeking a solution to (41) for any function f in some space is equivalent to 
seeking the inverse of the operator L in tliat space. If the null space 
of I. is sero (i.e. the solution of (45) is only the trivial solution) 
then I, is one to one and its Inverse may be defined. If (4b) has non- 
trivial solutions, L is not one to one and h ^ may bo defined, if at all, 
not uniquely on the range of h. The consistency vondltlon guarantees that 
f has no component in the null space of L, hence (with a little more 
work) that it is in the range of L. 

If we approach the determination of the Green's function as we did for 
the beam with simply supported endpoints then the attempt to solve the. BVP in 
(46) results In on immediate stumbling block, This system has a solution only 
if 


(a(P-Q), v(Q)) » 0 (50) 

for every solution v(Q) of (48). But, the delta function, in general, 
has components In the space, of solutions to the homogeneous problem, so 
instead we seek the modified Green's function which satisfies 

L g(P,Q) - <S(P-Q) - I u^(P) u^(Q) (51) 

i 

where the u^(P) are the normalised, non-trivial solutions of (43). 

We have subtracted the offending components of 6(P-Q) which lie in 
the nullspace of L. A solution to this system does exist. 


t3 
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If the Green's function Is known, the solution (47) of the boundary 
value problem is known. If b(F,Q) is not known (47) is an expression of 
the solution in terms of a compact operator which incorporates the boundary 
conditions. Approximate solutions may be computed from the eigenfunction 
expansions 

8(1’, Q) * I f' 

k«l ^k ^ ^ 


and 

0 to 

Kf - fjj g(P,Q) f((J) dQ » I <jSj^(P) 
* k*l k 


(53) 


where are the non-zero eigenvalues and are the corresponding normalized 
eigenfunctions of (41-42), which satisfy * Xj^ s{ij^ and B^(^j^) *• 0, 

1 < i < N. 


Substitution of expression (52) for f In (53) yields the following 
relation which will also he useful: 

00 

SCP^.Q) tUO.lV UQ. I rj2-*;,(V • 

•' k»l k 




(54) 


Alternatively, other numerical methods may be applied with greater 
convenlenee to (47) then to the original boundary value problem (41-42) , 
because of the superior properties of Integral over differential operators. 


5,3 The General Shape Control Problem 

In this section we define a general control problem and a general 
estimation problem corresponding to large space structures. We then solve 
these problems using the results of the last section. In the next section 
wo will give specific examples o£ problems and their solutions. 

The control/sensor mechanisms for large space structures will probably 
be located at discrete points P^, 1 £ i £ m, along the structure, rather 
than continuously. Thus the general dynamical model for the control 
problem is 

m 

bu » J f, 6(P-P.) (55) 

i»l ^ 

B^(u) - 0 , 1 < j 1 N 
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(56) 


whi‘re uCP) is the nhape, h 1b a linear differential operator as before, 
iB a force to be applied at the position P^, and (56) denotes an appropriate 
set of boundary conditions. 


l.et ^(i)) be the desired Bhape of the space structure, and define the 
quadratic criterion 


m 


■)(F.U) . I f/ + A 


(^(Q) - u(Q))^ dq 


(57) 


as a measure of performance, 
p • (fj . f . 


The constants q^ are arbitrary weights and 


The control problem is to determine the vector of forces F* whicli 
together with the corresponding solution u* of (55-56) minimises .1 over 
all admissible sets (F,u). 


Tlie solution of (55-36) is given by 


u(P) « 


m 

bO’,1!) 1 1 t. d(! 


l»l 


m 


I f. 8(F,F.) 


(58) 


l««l 


where g(l’,Q) satisfies (46), Substitution of (58) into the criterion (57) 
yields the criterion 


J(P) 


-• !U yx if 0 

i ^ 'll Jq ^ 

i*l i*l 


The constrained optimization problem (55-57) has become, the simpler problem 
of minimizing a function of m unknown constants without constraints. 

Solving simultaneously the equations 

A 

14'’* * 0 1 < i < m , (60) 

di^ - 


we arc lead to the following necessary condition for an optimal solution 

F* « (e,* ... f *)^; 

'1 m ' 

(Q + A) F* « 8 (61) 
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The m X m matrices Q and A have coefficients 


«i-J) 



jR 


B(P^,Q) g(P^,Q) dQ 


( 62 ) 


and the m dimensional vector B has coefficients 

B » [ g(P.,Q) «i»(Q) dQ . (63) 

Once the optimal forces are determined, the optimal shapa u’*' is given by 
(58). 

It is interesting to note at this point that the necessary indications 
for the shape control problem require only thct Q be positive semidefinite, 
since A is positive definite. That is to say, in contrast with the full 
state optimal controllers in the time domain, placing a zero weighting 
on a particular control will not result in unbounded control forces. A 
simple example will serve as the explanation of this phenomena. 

Consider an elastic beam pinned at both ends with a control force 
locateu it the center. If it is desired to bend the beam Into any shape 
whic)> is not symmetric with respect to the applied force, no amount of 
control force will accomplish the task exactly. That is to say, even with 
a zero weighting on the coat of control, a bounded value of control force 
will come closest to producing the desired shape. 


5.4 The General Shape Estimation Problem 

For the estimation problem we assume the shape u(P) satisfies the 
boundary value problem 

Lu » f , B^(u) = 0 , 1 i. 3 1 N , (64) 

where f(P) is an unknown function representing disturbances or errors in the 
model. Sensors placed at the positions P^, 1 £ i £ m, yield the 

observations 

y^ = u(P^) + v^ (65) 


43 


where v. is an unknown constant representing inaccuracy in the observation 
i 

at P^. Let V « (v, ... V ), We define the performance criterion 
i 1 m 


1 m n 1 

J(V.£) - 2 I 

i-1 


m 


rCQ) dQ 


( 66 ) 


2 I ‘■H - 'll ■*'2 

i»l 




£ (Q) dQ . 


The estimation problem is to determine the pair (u*,f*) which jointly 
satisfy ''64-65) and minimize the criterion (66) over all admissible pairs 
(u,f). 


The solution to (63) is given by 


u(P) 


g(P,q) f(Q) dQ 


(67) 


where g(P,Q) again satisfies (46). Thus 


\ = ^i - 




g(P^,Q) fCQ) dQ 


( 68 ) 


We substitute (68) into the criterion (66) , v/hich produces the criterion 

A 1 ^ * 9 i f 9 

J(f) “I I g(P^,Q)f(Q)dQ)^ ^i n ^ 

i=l 

A 

The problem is now to minimize the functional J without constraints, A 

A 

necessary condition for a minimum of J at f* is that the differential 
3J(f*,h) - 0 - J q^(y^ - L g(P^,(Q)fA(Q)dQ)(- ^ g(P^,Q)h(Q)dq) 


i=l 

r 


f*(Q)h(Q) dQ 


for all admissible variation h. Thus it may be concluded that 
m 


f*(P) = I q^ g(P,P^)(y^ - u*(P^)) 
i=l 


(70) 


(71) 


Substitution of this relation into (25) yields the optimal shape estimate 


m 


u*(p) = I [q^(y^ - u*(P^)) 


i=l 


g(P.Q) g(P.,Q)dQ] . 


(72) 


Note that u*(x) is expressed in terms of the unknown discrete shape 
estimates u*(P^). Let 


X = (u*(P.) . . . u* (P )) 
i m 


( 73 ) 
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and 

Y - (yj^ ... y„)''' (?•') 

Evaluation of (72) at x » Xj, j - 1, m yields the fallowing necessary 

condition for the vector x: 

(I + AQ) X » AQY . (75) 

where A and Q are the matrices of coefficients (62). 

Once the vector X has been determined the optimal shape estimate is 
given by (72). 

In the case the Green's function is not precisely known, or is 
unwieldy, approximate solutions may be obtained using the eigenfunction 
expansion (54) for the coefficients 

5.5 Examples of Static Shape Estimation and Control 
Case I Shape Control for a Simply Supported Beam 


Consider the problem of controlling the static deflection of an elastic 
beam of length i (fig. 12). Define, a coordinate system such that the x-axis 
passes through the endpoints of the beam, with one end at the origin and the 
other- at X - Suppose control is to be Implemented by means of transverse 


forces at positions x^, i 5. i £ m, where 0 < Xj^ < x., 


. . < X < i. 

m 


At each point x denote the deflection by u(x) . Assuming no net 

tensile force on a cross-section, the shape of the beam is governed by 
the differential equation 


i f 6(x-x.) 
i=l 


(76) 


The ends of the beam satisfy the boundary conditions 
u(0) = u"(0) = 0 u(i) ■= u"(i) = 0 


(77) 


Let 'jj(x) be the desired shape of the beam. As a measure of performance we 
define the criterion 




( 78 ) 


T 

where F ’s the vector of forces (f, ... f ) and q. are non-negative 

1 m ‘t 

constant weights whose values are optional. 
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The object is to dutetmine the set of forces wliich together with 

the solution u’*‘(x) of (76) minitnizes (78) over al3 possible pairs (u,F). 


The solution of (76) is given by 


m 


u(x) - I g(x,x,)f. 
i»l 


(79) 


where g(x,?) is the Green’s function associated with (76,77) which satisfies 


dx^ 

g(0,5) - g"(0,5) - 0 g(i,F,) - g"(S.,f,) - 0 . 


(80) 

(81) 


The Green’s function represents the natural response of the beam to a unit 
impulsive force at X » The existence and uniqueness of the solution 
to (80,81) follows from the fact that the associated homogeneou.-i system 

,4 

« 0 v(0) « v”(0) « 0 v(£) » v"(£) * 0 (82) 

dx“ 


has the only trivial solution. The solution of (80,81) is 


1 

; ■ 6£ 


0 K 


(x''-2£x+5^) 

Jr V 0 

6S. 

^ A A* 




(83) 


Figure 13 dis\>lays the Green’s function which corresponds to impulsive forces 

S 

at positions 5 * n (-g-) , n * 1, 7. 


The solution of the control problem follows, 
solution (79) into the criterion (78) yields 
A 1 m 


Substitution of the 


J(F) 


1 


i“l 




m 

n 


1*1 


g(x,x^)f^-\jl(x))' 


dx 


(84) 


The problem of minimizing tl\e criterion (78) subject to the constraints 
(76-77) has become the. problem of minimizing a function of m unknown 
constants without constraints, A necessary condition for J to liave a 
minimum at F* is 

A 

~ “0 1 < i < m (85) 

C j 
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ThiH condition bccnmoB 
m ei 

0 R(^»^i)«Cx,Kj^)dx " i{i(x)R(x,x.p dx 


( 86 ) 


ir we define 
tl 


and 


l\ »* . 

1.1 j 


o 


i»(x,x^)r,(x,x^) dx, 1 < i, J < in , 


rk’ 


b «* 
"i 


I|j(x) |*(x,x.p dx, 


1 < 1 m , 


(67) 

( 88 ) 


then the neceBsnry condition fod' a uiininmm of J at F* is that F* BatiHfy 


(I + AQ)F « A q B 

where Q ie the m x in diagonal matrix 


(89) 


Ah 


[ 


\ 


(90) 


q ? 
‘ill 


A is the m x in matrix with coeffleienta (87), and B is the lu dimensional 
vector with coefficients (88). 

the s imply supported beam 

(1) Compute the constants a^^ and defined by (87-88). Define Q,A,B. 


m 


(2) Solve (89) to get F*. 

(3) The optimal shape n*(x) = ^ f,*g(x,x.). 

i=l ^ ^ 


2tt k 

Figure (14) displays the optimal shape vs. the, desired shape ijj(x) » sin 
the second mode of the system (76,77), for two actuators at 1/42. and 3/42. 


Case II . The Control Pi'oblem for the Pinned-Free Beam 

A modification of the control algorithm is necessary if the system 
has rigid body modes, as is the case with the pinned-froe beam., 

The beam with one pinned and one free end point satisfies the diffcreiitial 
equation (76) with boundary conditions 

u(0) = u”(Q) * 0 u"(£) « u'" (2) « 0. (91) 
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Wq will ngaln use cho pcrforaonca crltotlon (78). TUc object is to 
determlse the set of foteea (f^) which tosetUet with the solution u(n) 
of (75, Si) minimices (76) over nil possible poles ((E^) , n) 

The system (76,91) has the rigi<l body mode == (normalized) 

Physically this means the beam can hove a non-zero slope or tilt as a rig c 
body. Mathematically it means that the corresponding homogeneous system 


^ = 0 v(0) = v"(0) ■= 0 v"(i) = v"' U) = 0 

dx'' 


(92) 


has the non-trivial solution Uj(x) . Thus the system (76,91) has a solution 
only if tl'G inner product 

m /T“‘ n - „ (93) 


Tw ^ ni 

(j^E,6(x-x,),u,) = J,Vi " “ • 


The additional constraint (93) must be added to the problem of determining 
the optimal control forces. 

A solution to (60) with plnned-Erce boundary condltloas does not exist 

I »• (xfv {'^ II ^ is not zero. The ''modified Giecn s 
because the inner product (<S(x-0, u.^ is not zero. 

function which is appropriate to cho system (76, 9i) satisfies 


= 6(x-t) - -i; 

rt J 




cl 

dx^ 


(9^0 

(95) 


we make the additional requirement that r„^(x,C) have no component in the 
subspace spanned by the rigid body modes. 


(Rr 


(x,e),u,) = dx = 0 


(96) 


The modified Green’s function which satisfies (9A-96) is given by 


, , , ,33a . r+x^. _ g„+x_,) 

gjxyf,) = (j^ + •- 


2. 2 
, +i 


3 


^0.?, 


1 3 3 

+ 0<x<5 

2 3 

4. 4— £!<X<R' • 

2 6 ^ 


(97) 
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Condition (96) guarantees that gj^^Cx,?) is synunetrlc and of minimum norm 
among all solutions of (94,95), 

The Green's function (97) represents the response of the pinned- free 

beam to one of a set of unit impulsive forces which satisfy (93), Figure 15 

displays the Green's function for impulsive forces at positions n (|) , 

8 

n 1,,,,,7, 


The solution of (76, 91, 93) is given by 


m 

■j(x) ° I f. g (x,K.) (98) 

1=1 

We solve (93) for f^^ in terms of the other forces and substitute that 
expression together with (96) into the criterion (78), which results in 

*>imo * 

J(F) = ^ I 


1=2 1 


1=2 


1 ( ^ "J X • _ 

2 j ^ ” x^ fi(x,x^))-ij;(x))'^ dx 


(99) 


where F is the vector (f,,.,,f ) 

* m 


T 


Again, the optimization problem is reduced to one of minimizing a 
function of unknown constants. 

The necessary condition for a minimum at F* is 


8J 

3f. 


(F*) 


2 < 1 < m , 


( 100 ) 


The conditions (100) result in the following algorithm. 


(1) Toinpiite the m dimensional vector B and m x m matrix A whose 
coordinates are 

r£ 


= 8,„(x,x^) i/((x) dx 

• rt 


*ij 


dx . 

in “J 


( 101 ) 

( 102 ) 
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(2) Computa the (m-1) dimensionaX vector B and (ro-1) x (ni-X) matrix A 
whOBC coordinates are 

A X, 


•^1 “ *^ 1+1 


Jdk b 

X, “l 


”ij * .. 2 




(103) 


^i+l,j-fl " ®l,i+l X 


fl+1 


H+1 


"ij+i x^ 

A 

Lot Q be tha (m-1) x (m-1) diagonal matrix 


^2 


r 


r' 




(3) The vector F* of optimal forces satisfies 

A A A A 

(Q + A)F* * B . 


(104) 


(105) 


(106) 


The optimal force is found from (93) . 

m 

(4) The optimal shape u’'f(x) » V f,*g (x,x.) . 

1«1 ^ ^ 

Since the optimal shape u* is a linear combination of (Jreen's functions 
which satisfy (96), it will have no component in the subspace of the rigid 
body mode. If the desired shape ^(x) does have sucli a component, that is if 
(t|;,Uj^) is not xero, the optimal shape will approximate the shape 

Tp(x) - (ili,Uj^) Uj^(x) . (107) 


That is, it will approximate the desired shape minus its component In the 
subspace spanned by Uj^(x). 

2 

As an example. Figure 16 displays tlie desired shape \p(x) - Hx“X , the 

3 2 

shape which approximates £x-x , and the optimal shape plus the missing 

rigid body mode component 3. £xo 

4 
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Those components of the doHircd shape In the subspace spanned by rigid 
body modes must be added by the attitude control system. A shape control 
system constrained to satisfy the boundary conditions cannot affect these 
components. 


Case IV. The Shape Estimation Problem 

To illustrate the shape estimation algorithm we consider a simply 
supported beam of length i and unknown shape u(k) , which satisfies 

,4 

— !:} * f(K) on 0 < s < £, (108) 

dx^ “ ” 

and 

u(0) * u"(0) * 0 u(jl) « u"(fc) » 0 

The function f(x) represents minor model luaccurncies or random disturbances 
acting on the beam. 


Assume sensors at positions 


0 < < . . . < produce 


observations 

\ ^ -1 ^ 1 • (109) 

As a measure of the accuracy of shape estimates we define the criterion 

J(f,u) » I (x^-u(x,^))^q^ 4* I p f^(x) dx . (110) 

i»*l M! 

The object is to determine the function f>^ which together with the solution 
u* of (108) minimizes (110) over all possible pairs (f,u). 

The solution of (108) is given by 

^jL 


u(x) 


g(x,c) f(?) d5 


( 111 ) 


•' 0 

where g(x,5) is the Green's function (83). We substitute (117) into the 
criterion (110) ; resulting in the criterion 


J(f) 


T m 

2 p " 
1»1 


g(K^»5)f(OdO^ +-| 




f(5) dC. 


( 112 ) 


The estimation problem has reduced to one of minimizlnR (112) without 
constraints. A necessary condition for J to have a minimum at f* is that 
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the Frechet d' ferential 


A | 1 } 

DJ(£,h) « I q.(K. “ g(x.,C)l^*(e) d£) 

i-l Jo ^ 

rl fi 

(‘ 

'0 •'0 

for all admissible variations h. This implies 

m 

f*(0 * I qig(x^»5)(8j^-u*(xp) . 


C“ 8(x.,5)h(g)d5) + f f*(F,)h(g)dC - 0 

^ o J n 


Then 


i»l 

m 




Lot 


and 


u>^(x) « I g(x,f,)g(x^,g) df,. 

i*»l J 0 

X « (u*(x^) ..« u’^(Xj^j))'^ 


(113) 

(114) 


Z * *** * 

Evaluation of (114) at x ■ x^ and regrouping of terms yield the following 
neeessory condition for the vector X; 


(I + Aq) X * AQE 


(115) 


whore A is the matrix of coefficients (87), Q is the diagonal matrix (90). 


The shape estimation algorit hm . 

(17 Compute the elements of the matrix A given by (87), and define 

X,Q,^» 

(2) Solve the system (115) for the vector X. 

(3) The optimal error estimates are given by (41) and - u’'f(x^), 

1 < i < m. 

(4) The optimal shape estimate la given by (114). 


This algorithm is equally valid for the static beam with other boundary 
conditions, provided the appropriate Green’s function is usecL 


Figure 17 displays the optimal shape estimate versus the actual shape 

. /TTXv , 1 j /2irXx 

sin (j“ ) + sin (~) » 

11 3 

for three exact observations at j it, and ^ 5,. 
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5.() Approx^^.inations for Shape Estimauion and Cmxtrol 


The following approximations arc purtieular^- 4 . nimble when the Hyatein 
models are more complicated so that the Green’s functions are difficult to 
compute analytically. They are based on the assumption that the (Jreen'vS 
function is symmetric and of minimum norm, whiel» fallows if the boundary 
value problem is self-adjoint. 


n m 

ii 




k«i a: 


CO 


i/<x)g(x,x^)dx - I "X 
k**l 


(lib) 


whore x^’s are the actuator (or sensor) positions, Xj^'s are the non-zero 
eigenvalues, and the are the corresponding normalized eigenfunctions 
of the associated boundary value problems. 


Simulations of approximations based on the first term 
expansions generated approximate optimal shapes which were 
indistinguishable from the optimal shapes, and numerleaily 
third or fourth significant figure. 


in each of these 
visually 
accurate to the 


I. 
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Section 


6„0 Experimental Veriricatlon of Distributed Control ConeeptH 
b,l introduction 

The previous sections of this document have contained new developments 
in the theory of the control of large space structures. In some cases, 
computer simulations of these new developments have been performed to verify 
the ivnalytical results. One vital phase of this research and development 
program still remains to be presented, experimental vorif iciition. 
l.aboratory veri I lent ions of the analysis and computer simulation are required 
to fill the gap prior to the final, flight project, stage of the large 
structure work. 

An experiment employing a pinned~free flexible beam has been constructed 
to demoustrate and verily several facets of the control of flexible structures. 
The desired features of the experiment are to demonstrate active shape control, 
active dynamic control, adaptive control, various control law design approaches, 
and associated hardware requirements and mechanixaUlon difficulties. This 
section contains the analytical work performed in support of the facility 
development, the final design specifications, control law synthesis, and 
some preliminary results. 

The flexible beam v»>as chosen for this experiment for being a simple, 
continuous structure with many of the dynamic characteristics that are 
representative of general large space structures, including, infinitely 
many vibratio.i modes, a rigid body mode, and many, "low" frequency vibration 
modes . 


The selection of the flexible beam alvSo resulted in some minor limitations. 
The flexible beam does not Jiave repeated eigenvalues, however, by orthogonaliily 
of the corresponding eigenvectors, these modes may be distinguished spatially, 
if not by frequency domain metlmds. Secondly, totally free boundary- 
conditions are not possible in a ground based experiment, IVo beam support 
configurations which replace the rigid body mode with a lov^.* frequency 
pendulum mode were considered as alternatives and are shown in figure 18. 

hO 



UJ 


(b) 


Figure 18, Beam wSupport Configurations 



Figure 19. Beam Suppox"t Structure 
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Figure 18a shows the beam supported by several cables In the horizontal 
direction and figure 18b shows the beam hanging vertically from a pinned 
support. Although the configuration in figure 18a more closely approximates 
the totally free boundary conditions, the complexity of distributing the 
weight of the beam uniformly and the difficulty with constraining the unwanted 
degrees of freedom in 18a resuJ ted in selection of pinned-free hanging beam 
of lb as the final configuration. 

One further compromise had to be made in the sensing/actuation phase 
of this experiment, A totally free space structure is constrained to sense 
with respect to (and react against) itself, or inertial space. Because the 
major objectives of this experiment are to demonstrate control technology, and 
not sensor/actuator technology, sensing and actuation of the beam are both 
performed with respect to an external frame. 

6.2 Dynamic Analysis of Flexible Beam Facility 

A schematic of the beam and its support structure (tower) as they 
are being erected is shown in figure 19. The tower is constructed of 
aluminum angles and is twenty feet tall, two feet deep in the stiff direction 
and one foot deep in the compliant direction. The weight of the tower is 
two hundred pounds. With the sensor/actuator mounting brackets (figure 20), 
beam, sensors, actuators, and electronics, the total weight is about three 
hundred pounds. 

Shake tests were performed on the tower to determine its resonances, and 
if they might interact with the control of the flexible beam. The results 
are given in oi*der of increasing frequency in Table II. 


Table II Tower Resonances 


Frequency 

O^z) 

Mode 

Direction 

6 

cantilever 

compliant 

10 

cantilever 

stiff 

27 

pinned-free 

compliant 

35 

pinned-free 

Stiff 

45 

free-free 

compliant 

63 

free-free 

stiff 
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For the final configuration, only the natural frequencies in the stiff 
direction are of interest because there is very little coupling between 
the perpendicular directions. 

A schematic of the sensor/actuator mounting bracket is shown in 
figure 20. Either the sensor, or actuator, or both may be mounted on a 
single bracket, and the brackets may be mounted at any of the stations 
located at six inch intervals along the beam. 


A dynamic model of the hanging, pinned free beam is required for 
control system design. Two modeling approaches will be examined here; 
first, an analytical approach, and, second, a finite element model. 


Temporarily setting aside the effects of the beam being in tension 
due to gravity, the partial differential equation of motion for the elastic 
beam with constant mass and stiffness per unit length and the appropriate 
boundary conditions are 


p ilx + 121 iZj » f 
Ot 3x 


(117) 


.2 

y(0,t) »• 0 — ^ (h»t) = 0 
Ox 


, 2 3 

(0,t) « 0 ^ (L,t) » 0 

Dx 3x 

It is straiglitforward to show by assuming an eigenvector decomposition 
of the solution that; the eigenvalues and normalized eigenvectors are 


+„(C) 




L 


0 ) » 0 
0 


(118) 


12 2 2 
cosh^X cos^X 

I \ . L n n 

"\l ;7' 2 ~ 

cosh X - cos X 


sxnnX C sinX ? 
n , n 

T 


cosh X 


n 


cosX 


n 


v^rhere 


w 


and 


and 


EIX 


pL 


tanhX “ tan X 
n n 


E - 
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Using nn asymptotic approximation (good for n 1) for (118) 




/? 


oj » 0 

o 


♦„(5) = )j; 


slnX g slnA K 

_ _ ^ n_ ^ n 

ciisTTr’ c’bvr*'’ 


(4n+l)n 


( 119 ) 


n * 1,2,3, , 


The dynamic equations for the modal amplitudes, become 


q + « q 
‘n n ‘n 




n » 0,1 ,2 . . . 


( 120 ) 


For the specific case that fOi,t) 
concentrated force applied at f, * 


F(t) 5(^^), i.e., a spatially 


( 121 ) 


A graph of the first four mode shapes is shown in figure 21. 

The rigid body motion of figure 21 is a zero frequency eigenvalue. Of 
course, for the plnned-free beam hanging under the influence of gravity, 
no such mode exists. Rather than being a zero frequency, rigid body mode, 
the actual dynamics is a low frequency, pendulum-like bohavlur. In fact, 
gravity interacts with all of the modes to some degree. To determine which 
behavior (tension or elastic) dominates the various modes, an Independent 
dynamic analysis of a hanging free string will bo performed. 


The partial differential equation of motion for a hanging string, 
with its boundary conditions is 


a^y 3 3y 

SS ».- g— CTJn .(Y V* mirntT 111,1 

3^2 ax 


y (0) finite 


( 122 ) 


y (h) *= 0 
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The solution to (122) is given by 

4, , j (u /H) 

^n 0 ' n g 


where 




/X) 


0 . 


(123) 


Again, an osyroptotic approximation to the eigenvalues can be obtainet) 
“n 2h 2 


If the true mode shapes of the exact system are obtained from 124 
2 4 

p ilZ + El 5-.-^ “ n g (t»“X) * f (124) 

Dt"^ 3x 

it can be shown that the squares of the exact eigenvalues are approximately 
the sum of the squares of the string and the beam frequencies with the largest 
error occurring when the two frequencies are equal . 

Using the experimentally obtained beam parameters appearing in Table m, 
vhe approximate modal frequencies and the experimentally determined modal 
frequencies can be compared. These results are contained in Table IV. 

Since the exact analytical solution to the differential equation of 
motion is tiot known, a discretized, finite element model is developed, Tlie 
beam is divided into N segments as shotm in figure 22. The displacement 
and slope at the ends of each segment arc specified, as coordinates. Mass 
and stiffness matrices can be defined for each element, and these assembled 
to create overall mass and stiffness matrices for the system. The differential 
equation of motion is thereby replaced by a matrix elgenproblem, and arbitrary 
accuracy can be obtained by considering smaller divisions of the beam. 

The stiffness matrix K for the finite element shown in figure 23 is 
defined by 

L! = l/2x\x ( 12 :) 

where U is the potential energy of the segment and 

X « (x^, X2» x^) 
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(126) 


g 



x*0 


f ..-f- ..-f -..-f *-h -f 


0, HX, % 




Figure 22, Fiiilte Element Beam Model 


l''iiilLu Hiemuia 




f), EX 


x“Xj^+(Sx 



Figure 23. A Single Finite Element 
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gives the nodal clisplacementa and alppeB, Analytically U la given by 

U - in 


V 

^ (Ely"^ + pgU - x)y'^) dx 


X, 


Substituting 



we obtain 


U • 

M 


2,+(Sz 


(y*'^ + Zy’^) dZ 


wheru 


z, * 




nnd 


6z - - 5x 


( 127 ) 


(128) 


(129) 

am 


The displacement y(z) of the finite e]emt»nt is chosen to minimize the 
potential energy U of the element. This is equivalent to asHviming that 
the shape of the element Is unaffected by dynamic loading. 

Equating a variation of U to zero and integrating by parts, we 
obtain 

y"" *“ (zy') 0 (131) 

whose solution is 


y(z) » UqC^qCz) + C‘ 3 ^‘}>j^( 2 ) +«24’3(x) 


(132) 


where the a^’s are constants and the are linearly independent 
solutions of (131) . Though the analytical solution to (131) is known in 
terms of Bessel functions of fractional order, It is more convenient to 
define the functions os infinite scries. Tin' ■‘b^ir.e used in this 
analysis is 



(133) 


^^(x) • 1 

•» (z - Zj) 


+ 


(z *• S 5 ^)^' 


+ 


• f • 


iw 


(Z - Zj^)^ 2j^(z - Zji^)^ 

. ,„. — + ._-.Y|. — + ... 

(z - Kj)^ (z - Zj^)’ 

SI »"■ » -w C. # # * 

31 1 51 


with eoef fleitxntB of (z - z^) (Jutermlned by the recuraion relaeion 


mi 


a 


n 


^l‘'rt+l 


n+3 


(a + 1) (u + 2)(n + 3) (n + 2) (n + 3) 


n- 1,2,3 ,... (134) 

T1 h 5 functioiiH art* expanded about z * Zj to enauro rapid convergence of 
the HcrieH. 

Hitnllar BorleB expansloiiH for the dei‘lvatlvea of ore eafilly derived. 
Of eourae, in practice the sorioB ore truncated after a finite number of 
torniH, once convergence has been determined. 

Kquation (132) can be written in matrix form as 

y(z) * (}i^a (135) 

where 

■(i})y(z), »}>3^(2)» <}>3(2)) (136) 



Furthermore, a can be related to q. Writing 


T 

n ** ^2> ‘la* ‘Ia^ 


I I 

* (y(xj^), y(i<jL)» y(j<x 

and substituting (135) into (137), we obtain 
q * Aa 


(137) 


(138) 
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where the matrix A is given by 


A >* 




ij)(Xj^ + fix) 
ij)^(Xj^ + fix)^ 


Equations (138) can be rewritten as 


( 139 ) 


a « Bq 


(140) 


witli 


B » A 


-1 


Substituting (140) into (135), and that in turn into (128), we obtain 


T, Tnl* 

U = "~2 ^ 


^z^ +fiz II llq, , , m 

' ((() (|) + Z(f> (}) ) dzBq . 


(141) 


Identifying terms with equation (125), we finally obtain an expression 
for the stiffness matrix: 


K = - pgB KB 
where 


K 


rz.+fiz II 11 ^ I ,m 

(tfi (|> + z<j> (j) ) dz 


(142) 


(143) 


The matrix K can be evaluated by Integrating (143) by parts twice and noting 
(131), which yields 

z.+fiz 

(144) 
"l 


<' 1 1 1 ni 111 fp 'p 

K“((|)(j) “(}) cp +Z{j)’(j)) 


The mass matrix M for this finite element is defined by 

rr 1 *T , • 

T = J q Mq 

T is the kinetic energy of the elements 


(145) 


For the beam element. 


T = ^ 
^ 2 


fzx+fiz .2 
' py^ dx 


^ 1 


(146) 
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An analysis analogous to that of the previous section yields 


M 

and 


rp^ 

Y 


M 


f^‘1 ^ n’ 

<j) (|) dz 


and where y» B and (p are defined as before. 


(U7) 


(148) 


The integrals in (148) are evaluated by expanding the integrands in 
series and integrating term by term. Three terms were taken, since it 
can be shown that truncation of the series after three terms results in 

A 

less than .2% error in the elements of M for a 20 division finite element 
model. 


Tlie beam shown in figure 22 is divided into N finite elements and is 
represented by 2N + 2 coordinates x^. The overall stiffness matrix K is 
defined by 

(149) 


1 T 

U « x\x 


Equating U to the sum of potential energies of the finite elements, and 
matching displacements at nodal points, we find that 

k 


K 


a. 

11 

k 

12 

13 

k 

*^14 

0 

0 

1 

12 

I. 1 

^22 

k 1 

23 

k 

*^24 

0 

0 

1 

13 

1 1 
^23 

I 1 

^33 ■^^11 

^4 '*■^12 

1 2 
^12 

k ^ 
^14 

1 

14 

k ^ 
24 

k ^+k ^ 

34 ^'^12 

^4 ■^‘^22^ 

2 

“'23 

• « 

k 2 

^24 


0 

k 2 

13 

k 2 
23 

k33^+Kj3 

2 3 

*^34 '*■^12 


0 

k ^ 
41 

k ^ 
*^42 

“ 34 ' 

1 ^ 44' «22 
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A till) till 

where K. , donates the (l,j) component of the ft finite clement 

Xj 

stiffness matrix. Similarly, the mass matrix for the beam is given by 


o. 

11 

4. 

"’12 

a. 

'"13 

a. 

"’14 

0 

0 


1 

12 

1 

‘"22 

1 

"’23 

1 

*"24 

0 

0 


1 

13 

1 

"’23 

1^ 2 
’"33 +’"11 

"’34^-^"'l2^ 

2 

"’13 

2 

"’12 


1 

14 

1 

"’24 

1 2 
’"34 ‘^"’12 

1_^ 2 
"’44 +>"22 

2 

"’23 

2 

"’24 



0 

2 

"’13 

2 

’"23 

2 

'"33 

3 2 

■^’11 ‘"34 

•*■"'12 


0 

2 

"'l4 

2 

"’24 

2 

’"34 

^3 2 

■^"’21 ‘"44 

*"22 


(151) 


The dynamic equation of motion in matrix form becomes, 


Mx + Kx ® Q 
where 








.) 


gives the forces and torques applied at noda] points. 


(152) 


Tlu> boundary condition y(0) = 0 is satisfied by specifying the 
thereby reducing the order of (151) by one. 


Equation (152) can be written in state-variable format as 


0 

I 1 

• 1 f 

1 ,3il + i 0 . 

-m"^k 

— , 

0 j 

1 

] 

' ' M‘^Q. 


(153) 


The normal mode shapes and frequencies of the free vibration problem are 
the eigenvectors and eigenvalues of the state-variable dynamics matrix. 


) 

0 

1 

-1 


-M \ 

0 


(154) 
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Finding the eigensystem of (154) is not the most economical method of 
finding the solutions to (152) , particularly since the matrices M and K 
are symmetric and banded; however, the solution to the closed-loop system 
once feedback control is incorporated is most easily solved in state 
variable form. 


A Fortran program was written to generate the finite element stiffness 
and mass matrices, assemble these into the overall stiffness and mass 
matrices, and solve the eigensystem of the state variable matrix (154), 

The normal mode frequencies for N=10 and 20 divisions of the beam are 
summarized in table IV, The approximate analytical results and experimental 
results are included in table IV for comparison. These numbers are based 
on the following determination of the beam properties. 


Table III Beam Characteristics 


Material 

Length 


Stainless Steel 
149.875 inches 


Width 
Thickness 
Linear density 
Stiffness 


6 inches 
1/32 Inches 
0.644 Ib/ft p 
424.352 lb-in 


Table IV Normal Mode Frequencies (Hz) of Beam 


n 

Analytical 

10 Divisions 

20 Divisions 

Experimental 

0 

.301 

.308 

.308 

-.34 

1 

.728 

.755 

.755 

.75 

2 

1.27 

1.38 

1.38 

1.37 

3 

1.98 

2.21 

2.21 

2.15 

4 

2.92 

3.25 

3.24 

3.16 

5 

4.08 

4.51 

4.47 

4.38 

6 

5.49 

6.00 

5.93 


7 

7.13 

7.76 

7.62 


8 

9.03 

9.79 

9.55 


9 

11.16 

12.04 

11.73 


10 

13.54 

15.92 

14.15 
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In the near future, the following Improvements to the model will be made: 


• The actuator linkage masses will be included. These masses 
add to the mass matrix at nodal points in a straightforward 
manner. 

• Damping will be modeled. A common problem with high performance 
control systems is Instability of lilgli frequency modes. Since 
damping is Instrumental in the oCabQlty of these modes, it 
should be included in the model during the design process . 


6.3 Control System Hardware and Software 

In addition to the flexible beam, a variety of other components are 
required to complete the flexible beam facility. In the order they will 
be discussed, sensors, actuators, a microcomputer, microcomputer interfaces, 
and control software development, constitute a part of the completed 
facility. 

Both an optical position sensor and an eddy current position sensor 
reached the final selection phase. Many other possibilities were eliminated 
by the requirement of minimum sensor Interaction with the beam dynamics. 

In order to minimize the effects of external disturbances, tlie developmental 
period, and the final cost, an eddy current sensor made by Kaman Science 
Corporation was selected. 

The final actuator selection was a brushless DC torque motor manu- 
factured by Aeroflex Laboratories, Inc. With a three inch moment arm, and 
the appropriate mechanical linkages, the actuator has the capability of 
applying five ounces of force to the beam for a maximum one amp input, 
according to the manufacturer’s specifications. 

The purpose of the microcomputer in the control loop is to sample 
the sensors, pass this sampled data through a digital filter, and send 
the filtered data to the actuators. Assuming a general format for the 



digital filter 


\+i - * ’'k \ 

“k \ 

2 

The amount of computation is roughly n +np+nm multiplications. For 

the state vector larger than either the number of sensors or actuators, the 

2 

amount of computation is governed by the n term. Currently available 
eight bit microprocessors are capable of performing one fixed point double 
precision software multiply in one millisecond, and one hardware floating 
point multiply in one hundred microseconds. Due to the limitation imposed 
by the sampling theorem (i.e. sample two or more times per cycle) the 
maximum number of controlled modes with software arithmetic is four, and 
the maximum number using hardware arithmetic is eight. 

The microcomputer chosen for the control function is the SYM-1 by 
Synertek System Corporation. It is based on the 6502 microprocessor and 
has provisions for 4 K of random access memory (RAM) and 6 K of read only 
memory (ROM). Additionally, a KIMSI Interface/Motherboard by Forethought 
Products has been added for interfacing directly to S-100 products, specifically 
the digital to analog converters, and the hardware arithmetic. 

The physical interface between the analog sensors and actuators, and 
the digital microprocessor is accomplished through D/A and A/D hardware 
and appropriate buffer circuitry. Sensor sampling is performed under 
computer control. The twelve bit analog-to-digital (A/D) conversion is a 
successive approximation technique performed in software with the use of a 
Vector Graphic Precision Analog Interface Board (PAIB) . Similarly, the 
twelve bit digital- to- analog (d/A) conversion is performed in hardware on 
the PAIB. 

Sensor buffer /amplifiers were used between the sensor output and the 
computer. The circuit has a high input impedance and eliminates high 
frequency noise and DC offsets. The scale factor from position to voltage 
at the A/D converter is five volts per inch. The sensor bandwidth is 
greater than thirty hertz. 
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The actuator driver presents a high input impedance to the microcomputer 
D/A converter and eliminates the high frequency D/A conversion noise and 
DC offsets from the actuator command. Current feedback is used to eliminate 
the inductive effects of the torquer. The scale factor from the force 
applied to the beam to the D/A voltage output is one ounce per volt. The 
actuator bandwidth is greater than thirty herta. 

Software has been developed to Implement the general digital filter 
discussed in the section on microcomputers. The entire, program resides 
in 2516 erasable programmable ROM by Texas Instruments, and is located on 
the computer. When the software is initiated, the computer samples the. 
sensors, updates the state estimate, and ovjtputs the control. Data for 
the program is loaded into RAM and consists of G, and K, and the 
dimensions of these matrices. The exact sample period, T, (in msec) is 
given by 

T * 0,983 n^ + 0,963 np + 0.963 nm + 0.258 n 0.154 m + 0.597 p 
+ 0,010 m^- + 0.725 . 

A listing of the assembly language software can be found in the Appendix. 
6.4 Control Law Design 

A variety of control laws may be implemented using the. general software 
discussed in the section on microcomputers. The particular approach used 
in the initial control law design (and the example presented here), is 
an implementation of a discrete. Kalman filter using a sixth order 
estimator. 

The controller is based on the discretized version of the following 
continuous system 

x = Fx + Gu + rw 

u = C X (156) 

z = H X 4- V 
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where 

T • • *1 

X " ^o* ^1* ^1* ^2* ^2 


p *» 


3.55 


0 

0 

i 

i 

! 0 
\ 
t 

i 

i 0 


0 

0 

0 

0 


0 0 

0 0 

0 -1 

21.60 0 

0 0 

0 0 


0 

0 

0 


0 

0 

0 ! 

! 

I 


68.72 0 ! 

L j 

With a single position sensor and a single force actuator at the free end 
H « [1 0 1 0 1 0] 

g'^’ b r'^’ « [0 0.18 0 0.18 0 0.18] 

This control gains, C, and the estimator gains K are determined by 
minimizing the appropriate performance index J. In the case of the control 
problem 

(x^ A X + B u) dt 

and for this particular example 
i 1 


A = 


0 0 
1 

0 


I 0 0 

For the estimation problem 


, B = 0.0025 


J = 


T - 1 T - 1 

(w Q w + R ‘*‘v) dt 
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and for this particular example 


Q •* E(w w^) « 10 ^ (m/sec^)^ sec 

R » E(v v^) » lo"^ (m)^ sec 

For a six state estimator, one control, and one sensor, 

T « 0.049978 sec. 

Using a discrete optimal system synthesis DOPTSYS^^^^ algorithm, the 
following results are obtained 


0.950 

0.047 

-0.045 

-0.000 

-0.044 

-0.000 

-0.251 

0.910 

-0.085 

-0.033 

-0.053 

-0.019 

-0.019 

-0.002 

0.953 

0.048 

-0.018 

-0.000 

-0.064 

-0.084 

-1.144 

0.939 

-0.044 

-0.019 

-0.011 

-0.002 

-0.011 

-0.000 

0.904 

0.048 

-0.059 

-0.082 

-0.068 

-0.032 

-0.374 

0.896 


0.0438 

0.0265 

0.0183 

K 

0.0165 

0.0103 

0.0110 

-5.3744** 

-9.4553 
^ -6.5416 
-3.7488 
-3.0469 
-2.2154 

and the equivalent frequency domain eigenvalues of the closed loop system 
are given by 

-0.85 + 2.07 j -0.38 + 4.66 j -0.21 + 8.29 j 

Cq = 0.38 q = 0,082 == 0«°25 
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Figure 24. Controlled and Uncontrolled Beam Responses 
Amplitude =1/2 Actual Displacement 
Time Scale = 1 sec/cm 





Figure 25. Unstable Control System 

Amplitude = 1/2 A'ttual Displacement 
Time Scale = 1 Sf;c/cm 
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These macriceB ate converted to the sixteen bit fixed point format 
used by Che microcomputer and are loaded into the controller as data. 

The results for this controller can be found in figure 24. Figure 24 
contains the uncontrolled and controlled responses to both an initial 
condition error on positions and to an impulse force input. 

Figure 24a shown the open loop response to an initial position error. 

The majority of the response is the lowest frequency, pendulum mode. It 
has an open loop damping ratio of -2%, 

Figure 24b shows the closed loop response to the same input as In 
figure 24a, using the controller described previously. Notice the much 
faster decay rate of the closed loop system. The corresponding control 
force applied to the beam is shown in figure 24c. 

It is now possible to partially verify the control analysis. The 
predicted closet, loop damping ratio for the low frequency mode is 0.38. 

Using the maximum overshoot of the closed loop position response for 
determining the closed loop damping ratio, t... closed loop damping ratio 
of 0.40 is experimentally obtained. 

Figures 24d, 24e, and 24f show analogous responses of the beam to an 
impulse force disturbance applied at the free end. 

One additional featutj! present in the closed loop responses in figure 24 
is the high frequency oscillation appearing at the tail of each plot. 

This phenomenon is due to control/observation spillover into the first 
unmodeled mode at two hertz. 

The deleterious effects of spillover in control system design can be 
demonstrated in the following example. By decreasing the cost of the 
squared control by a factor of eight, a second control system can be 
designed and implemented. The results of this control system are shown 
in figure 25. 

In this case the control spillover into the first unmodeled mode 
is sufficient to drive the system unstable. 
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Section 7 


7.0 Future Work 

Future work in the distributed control of lurge space structures will 
be concentrated In two areas: (1) extending the distributed control 
analytical techniques to include general (FK) models and (2) expanding 
the experimental test facility. 

Distributed control analytical techniques will be extended by replacing 
the continuum models (if needed) used for understanding the theory, with 
more realistic FE models. Since the theory developed using PDE models 
will apply to any linear operator, the exte.nslon of the. finite elenient 
models will bo straightforward. 

The flexible beam test facility will be expanded tc an interactive 
facility. This will allow a general user to work at the computer 
terminal, design shape and/or active control systems, to Implement these 
control systems with the, microprocessor control system, and to obtain 
graphic outputs and chart recordings of the final results. 
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